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Abstract 

We derive a holographic dual for a gauged matrix model in general dimensions from a first-principle 
construction. The dual theory is shown to be a closed string field theory which includes a compact two- 
form gauge field coupled with closed strings in one higher dimensional space. Possible phases of the matrix 
model are discussed in the holographic description. Besides the confinement phase and the IR free decon- 
finement phase, there can be two different classes of critical states. The first class describes holographic 
critical states where strings are deconfined in the bulk. The second class describes non-holographic critical 
states where strings are confined due to proliferation of topological defects for the two-form gauge field. 
This implies that the critical states of the matrix model which admit holographic descriptions with decon- 
fined string in the bulk form novel universality classes with non-trivial quantum orders which make the 
holographic critical states qualitatively distinct from the non-holographic critical states. The signatures of 
the non-trivial quantum orders in the holographic states are discussed. Finally, we discuss a possibility that 
open strings emerge as fractionalized excitations of closed strings along with an emergent one-form gauge 
field in the bulk. 
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I. INTRODUCTION 



Extracting dynamical information on strongly interacting critical states of matter is in general 
a hard problem in theoretical physics. Fortunately, there are classes of strongly coupled quantum 
field theories whose non-perturbative dynamics can be accessed through dual descriptions which 
become weakly coupled when the number of degrees of freedom is large. 

One such dual description that has been extensively studied in condensed matter physics is 
the so-called slave-particle formulation yj-lj]. In this theory, a gauge redundancy is introduced in 
order to take into account dynamical constraints imposed by strong interactions. Unphysical states 
introduced in the redundant description is projected out by a dynamical gauge field. In the large 
N limit, where iV is the number of flavor degrees of freedom, the dynamical gauge field becomes 
weakly coupled and emerges as a low energy collective excitation of the system. 

The slave-particle theory may be viewed as a mere change of variables which allows one to 
compute dynamical properties conveniently, which could have been computed using a different set 
of variables albeit more complicated. However, the real power of the mathematical reformulation 
lies in the fact that it allows one to classify various novel phases of matter beyond the symmetry 
breaking scheme[4]. In particular, those phases that support emergent gauge boson possess subtle 
quantum orders that make them qualitatively distinct from the conventional phases. Because of 
the non-trivial quantum orders, the phases with an emergent (deconfined) gauge boson can not 
be smoothly connected to the conventional phases. Signatures of the non-trivial quantum order 
include fractionalized excitations and protected gapless excitations (or ground state degeneracy on 
a space with a non-trivial topology). 

The gauge-string correspondence is another type of duality QHZO. According to the duality, a 
class of D-dimensional quantum field theories is dual to a (D + 1) -dimensional string theory. 
The question we would like to address in this paper is : Do those phases that admit holographic 
descriptions in one higher dimensional space possess non-trivial quantum orders ? If so, what we 
call holographic states that can be described in one higher dimensional space can not be smoothly 
connected to the conventional non-holographic states. We claim that the answer to this question 
is 'yes'. The signatures of the non-trivial quantum order in holographic phases are the emergent 
space with an extra dimension, deconfined strings and the existence of an operator whose scaling 
dimension is protected from acquiring a large quantum correction at strong coupling in the large 
iV limit, even though the operator is not protected by any microscopic symmetry of the model. 
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The paper is organized in the following way. In Sec. II, we start by reviewing the slave-particle 
theory with an emphasis on quantum order in fractionalized phases. In Sec. Ill, we introduce 
a gauged matrix model which will be the focus of the rest of the paper. The model is general 
enough to include the U (N) gauge theory. In Sec. IV, through a first-principle derivation, we 
show that the matrix model in general dimensions is holographically dual to a closed string field 
theory in one-higher dimensional spaces. In Sec. V, it is shown that the partition function of the 
original matrix model can be interpreted as a transition amplitude between quantum many-loop 
states in the holographic description. In Sec. VI, we show that the holographic description has 
a gauge redundancy, and strings are coupled with a compact two-form gauge field in the bulk. 
Because of the compact nature of the two-form gauge field, topological defects for the two-form 
gauge field are allowed. In Sec. VII, we discuss possible states of the matrix model. Different 
states are characterized by different dynamics of topological defects in the bulk. If topological 
defects are gapped, strings are deconfined in the bulk, and the holographic state is stable. On the 
other hand, if topological defects are condensed, strings are confined, and the bulk description 
is not useful anymore. Suppressed topological defect in the holographic phase is responsible for 
a non-trivial quantum order which protects the scaling dimension of the phase mode of Wilson 
loop operators from acquiring a large quantum correction at strong coupling in the large N limit. 
We discuss the differences between the holographic and non-holographic states. The holographic 
critical phases can be divided further into two different classes. In the first case, there exist only 
closed strings in the bulk. In the second case, there are both closed and open strings, where open 
strings emerge as fractionalized collective excitations of closed strings. The latter state has a yet 
another quantum order which supports an emergent one-form gauge field in the bulk. Finally, 
we close with speculative discussions on a possible phase diagram, a world sheet description of 
deconfined strings, and a continuum limit. 

The present construction is beyond the level of identifying the equations of motion in the bulk 
with the beta function of the boundary theory. We construct a full quantum theory of string in 
the bulk that is dual to the boundary theory. The construction of the dual theory makes use of the 



fact that loop variab 
of matrix modelsj8- 



es associated with Wilson loops become classical objects in the planar limit 



the earlier works [12, 



The current construction of the string field theory is directly based on 
1311 . Compared to the the previous work on the U (N) gauge theory 111 311 . the 



present construction has two major improvements. First, the extra dimension generated out of the 
renormalization group flow is continuous, while the earlier construction produces a discrete extra 
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dimension. The infinitesimally small parameter associated with a continuously increasing length 
scale allows one to write the bulk action in a compact form in this formalism. As a result, one 
can readily take a continuum limit starting from a boundary theory defined on a lattice. Second, 
the earlier construction involves infinitely many loop fields in the bulk associated with multi-trace 
operators, which makes the theory highly redundant. In the present construction, the relation 
between single-trace operators and multi-trace operators are explicitly implemented. As a result, 
the dual theory can be written only in terms of the loop fields for single-trace operators. Because 
of these improvements, the dual theory takes a much simpler form, and this transparency allows 
one to uncover deeper structures in the theory. 



There a 



so exist alternative approaches to derive holographic duals for general quantum field 



theories [14J-120D. All these constructions including the present one are based on the notion that the 



extra dimension in the holographic description is related to the length scale in the renormalization 



group flow 112 1M25I1. 



H. QUANTUM ORDER IN FRACTIONALIZED PHASE 

In this section we review some of the key features of the slave-particle theory [jj-|3|] using a 
pedagogical model introduced in Ref. |2g]. We consider a model defined on the four-dimensional 
Euclidean hypercubic lattice, 

s = -*EE cos (^-^?) 

<i,j> a,b 

-ir^^cos^f + ^ + C). (1) 

i a,b,c 

Here 9°~ b, s describe phase fluctuations of boson fields defined at site i. Each boson carries one 
flavor index a and one anti-flavor index b with a,b — 1,2, N. < i,j > represents nearest 
neighbor bonds of the lattice. We assume that the phases satisfy the constraints 6 ab = —9 ba [5\]. 
With the constraints, there are N(N — l)/2 independent boson fields per site. The theory has 
global symmetry under which the boson fields transform as Of* — > Qf 1 + Lp a — ip b . 
In the weak coupling limit (K « 1), the model describes weakly coupled bosons. As the 
strength of the kinetic term t is increased, there is a phase transition from the disordered phase to 
the bose condensed phase. In the disordered phase, all excitations are gapped. In the condensed 
phase, there are (N — 1) Goldstone modes. (At the special point of K = 0, there are N(N — l)/2 
Goldstone modes due to the enhanced symmetry). 



In the strong coupling limit (K » 1), the large potential energy imposes an additional set of 
dynamical constraints, Of 3 + 9\ c + Of 1 = which is solved by a decomposition, 



e 



ab 



(2) 



Here 0"'s are boson fields which parameterize the low energy manifold. Note that these fields 
carry only one flavor quantum number contrary to the original boson fields. The new bosons are 
called slave-particles (or partons). The low energy effective action for the slave-particles becomes 



S 



(3) 



Note that this theory has a U(l) gauge symmetry, 



(4) 



This is due to the U(l) redundancy introduced in the decomposition in Eq. ©. Because of the 
gauge symmetry, the slave-particles can not hop by themselves. However, these particles can move 
in space by exchanging their positions with other particles. For example, in Eq. ©, the particle 
with flavor a can hop from site j to i as the particle with flavor b hops from itoj. In this sense, they 
can move only through the help of other slave-particles. One can introduce a collective hopping 
field Xij = Ylb e~ i ^ >i ~^ to characterize the amplitude of this mutual hopping. If we use this 
collective field, Eq. © can be written as 

S = -t Xi^~^. (5) 

<i,j>,a 

The magnitude of the collective field characterizes the strength of hopping, and the phase plays 
the role of the U(l) gauge field to which the slave-particles are coupled electrically. This map- 
ping from Eq. © to the U(l) gauge theory can be made more rigorous, by using the Hubbard- 
Stratonovich transformation[26]. Although the gauge field does not have the usual Maxwell's 
term, the kinetic energy is generated once high energy modes of the boson fields are integrated out, 
which renormalizes the gauge coupling from infinity to g 2 ~ 1/N. It is clear that slave-particles 
can propagate coherently in space only when the hopping field is 'condensed', and provides a 
smooth background. Since the hopping field is not a gauge invariant quantity, we need to be care- 
ful when we say that the hopping field is condensed. This notion can be sharply characterized by 
examining dynamics of topological defect. 



Because the U(l) gauge field is compact, monopole is allowed as a topological defect in the 
theory. The mass of monopole is 0(N) for a large N. Whether the slave-particles arise as low 
energy excitations of the theory depends on the dynamics of monopole. One can consider the 
following three different phases. 

1 . Confining phase 

For a small iV and small t, monopoles are light, and slave-particles are heavy. If monopoles 
are condensed, strong fluctuations of the phase mode of the hopping field confine the slave- 
particles. Only gauge neutral composite particles, which are nothing but the original bosons 
in Eq. (QQ), appear as low energy excitations. In this phase, all excitations are gapped. This 
phase is adiabatically connected to the disordered phase in the weak coupling limit. 

2. Higgs phase 

This is the phase which is electromagnetically dual to the confining phase. The slave- 
particles are condensed when t is large. As a result of the condensation of charged fields, 
monopoles and anti-monopoles are connected by vortex lines which produce a linearly in- 
creasing potential : monopoles and anti-monopoles are confined. One slave-particle is eaten 
by the massive U(l) gauge boson, and (N — 1) gapless bosons are left. These modes are 
the Goldstone modes. This phase is smoothly connected to the bose condensed phase in the 
weak coupling limit. 

3. Fractionalized (Coulomb) phase 

For a large N, the mass of monopole is large. When both the slave-particles and monopoles 
are gapped, the Coulomb phase is realized. In this phase, slave-particles are deconfined, 
and arise as (gapped) excitations of the system. They are fractionalized modes because 
they carry only half the flavor quantum number of the original bosons. Moreover, the U(l) 
gauge field arises as a gapless excitation. It is noted that the gapless excitation in this phase 
is not a Goldstone mode. It is not protected by any microscopic symmetry. Saying that 
there is a gapless gauge boson in a gauge theory may sound trivial. However, we have to 
remember that the gauge boson is nothing but a collective excitation of the original boson 
fields. The existence of a collective excitation which remains gapless without a fine tuning 
is actually something remarkable : someone who does not use the language of gauge theory 
would find the origin of the gapless collective excitation mysterious. It turns out that the 
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gapless mode is protected by a subtle order which is not characterized by any symmetry 
breaking scheme. This order, dubbed as quantum order|4(], is associated with suppression 
of topological excitation, monopole in the long distance limit. Formally, this order can 
be expressed as the emergence of the Bianchi identity dF = in the long distance limit, 
where F is the field strength for the emergent gauge field. The key features of the non-trivial 
quantum order is the presence of the fractionalized excitations and the emergent gauge field. 
Note that slave-particles are not gauge invariant objects. However, a 's become 'classical' 
in the large iV limit where non-perturbative fluctuations of the hopping field are suppressed. 
In this regard, fractionalization is associated with the emergence of an 'internal' space. 



TABLE I: 





slave-particle 


monopole 


low energy excitations 


Confining phase 


confined 


condensed 


gab 


Coulomb phase 


deconfined 


gapped 


cj) a , monopole, gauge boson 


Higgs phase 


condensed 


confined 


Goldstone bosons 



Table. I summarizes the physics in each phase of the boson model. Now, we switch gear to 
discuss about a matrix model and its possible phases. We will draw a close analogy between the 
quantum order present in the Coulomb phase of the boson model and a quantum order present in 
the holographic phase of the matrix model. We will see that the holographic phase has a distinct 
quantum order associated with the emergence of an 'external' space. 

III. MATRIX MODEL 

We start with a matrix model defined on the D-dimensional Euclidean hypercubic lattice, 

Z = [ dU e~ s[u] (6) 



with the action, 

<i,j> L 



S[U] = NM 2 HUtjUij) + N 2 V 



—W c 

N C 



(7) 
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Here i, j are site indices in the lattice with lattice spacing a, and Uij = is N x N complex 
matrix defined on the nearest neighbor bond < i,j >. Wc is Wilson line defined on the closed 
oriented loop C, 



W c = tr 



II u a 

.<ij>ec 



(8) 

where the product is ordered along the path. V[Wc/N] is a function of Wilson loop operators, 



V = -2>-" J2 J {Ci,~,c n }l[Wc k (9) 

n=l {Ci,..,C„} fc=l 

which is, in general, non-linear in the presence of multi-trace operators. Here J{Ci,..,c*„}' s arQ 
loop dependent coupling constants. This theory is invariant under the U (N) gauge transformation 
: Uij — > V^UijVj. Eq. © may be viewed as the partition function for a (D — 1) -dimensional 
quantum matrix model in the imaginary time formalism. 

To see that this model includes the usual U (N) gauge theory, we consider the following quartic 
action in Eq. © as an example, 

N2y = E {-NMl^ulp^ + Nv^ul^ulu^+v {tv(r u r ;j )}' 

<i,3> L 

-NJ^Wa, (10) 



where □ represents unit plaquettes on the lattice. Here Mq > 0, v > 0, v ' > 0, We assume 
that v is sufficiently large compared to J. The relative magnitude of M and M determines the 
shape of the potential for the matrix field. For small M , Uy = is the minimum, and the system 
is fully gapped. For large M , the low energy manifold is spanned by the matrices that satisfy 
Uijlljj = ul with u ~ ^(v~W ) " ^ n tn ^ s case ' tne ^ ow ener §y effective theory becomes the U(N) 
lattice gauge theory with the 't Hooft coupling A ~ ( Jtt 4 ) -1 . This theory can be viewed as a 'linear 
sigma model' for the U (N) gauge theory. Presumably, the gapped phase in the small M limit is 
smoothly connected to the confinement phase of the gauge theory. As M is increased further, 
the system can go through a phase transition to the deconfinement phase at a critical coupling M§, 
depending on the dimension. If the phase transition is continuous, we can take the continuum limit 
by taking a — > and M — > Mq such that the confining scale is fixed. 
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IV. GENERAL CONSTRUCTION 



In this section, we construct a holographic dual for the matrix model in Eg. © with general 
potential V in general dimensions. We will follow the idea introduced in Ref. il2ll where coupling 
constants are lifted to dynamical fields in the bulk space where the extra dimension corresponds 
to the length scale of the renormalization group flow. In the presence of multi-trace operators, 



this formalism becomes rather complicated[13] because one has to introduce independent fields 
for infinitely many multi-trace operators that are generated along the renormalization group flow. 
This issue is present even though multi-trace operators are not turned on initially, because they are 
generated at low energy scales in any case. To avoid this complication, here we express multi-trace 
operators in terms of single-trace one, by introducing a complex auxiliary field 4> c for each loop 
C (see Appendix A), 

Z = J dUd^dcpP* e~ s \ (11) 

where 

S 1 = NMMujjUij) + iV 2 4 0) (4 0) * - W c /N) + iVV[4 0) *]. (12) 



Here we dropped a multiplicative numerical factor in the partition function, which is not important. 
It is noted that Z is well defined although Si is not bounded from below as a function of 4>q and 
(f) c . This is because Si is complex and contributions from large negative Si is canceled because 
of rapid oscillation in phase. The repeated indices ij and C are understood to be summed over 
nearest neighbor links and closed loops, respectively. 



To perform a real space renormalization group II 12L 1271 28], an auxiliary matrix field is 
introduced in each link, 

Z = (ArV^)^ J d^d^dudtj e -& ( i3) 

where 

S 2 = Suy^A^-N^Wc 

+NM 2 tr(U} j U ij ) + NuMUijUij)- ( 14 ) 
Here Mi is the number of links in the lattice, and 

Suvi^AP] = N>hP^ + V[4>}\ (15) 



is an action for • We change the variables as 



U, 



-adz , 



Uij -\- Uij ) , 



Uij = e- adz (A Uij + Bu tj ), 
where a is a positive constant, dz is an infinitesimally small parameter, and 



A 



rn 



with 



M 2 

, B 



M 2 



m 



o2adz 



In terms of the new variables, the partition function becomes 



where 



Z = (N^mf 2 ^ J d^d^dudu e" 53 , 



S3 = S uv [^)\^)}-N<pfw' c 



(16) 



(17) 



(18) 



(19) 



+N 



(20) 



Here W c = tr 



and <f$ = e adzL °(j)p , where L c is the length of the loop 
C. The field with the large mass m has taken away a small amount of quantum fluctuations 

(q\i 

from the original field U^, which leaves an action for with smaller couplings (jr c . Therefore, 
we can interpret u^'s as low energy fields and u^-'s as high energy fields. 

Fluctuations of Uij renormalize the (dynamical) couplings for the low energy field u^. Inte- 
grating over u^, we obtain 



du e 



S 4 



(21) 



to the linear order in dz, where 



+N Mhx{u\ jUij ) 



N 



(22) 



n<ij> 6 c Ui i\ • m me tn ^ an( ^ me f our th terms, ij runs over all nearest neighbor 
links, and C\, C 2 are understood to run over all possible loops including null loops with the 



with w c = tr 
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FIG. 1: A loop with a self -retracting link splits into two loops (a), becomes shorter (b), or disappears (c), 
as the matrix field tty on the link is integrated out. In (d), two loops which share a link merge into one. 

convention 00. = 0, 00 = 1 and = 1 for null loops, where 0, refers to the null loop at 
site i. Here we regard null loops at different sites as different loops. By this, we can keep the 
combinatorics simpler. In the third term, F i3 -,[Ci, C 2 ] is a form factor that tells whether or not 
two loops C\ and C 2 are 'nearest neighbors' : Fy[Ci, C 2 ] — 1 if C\ and C 2 can be merged into 
one loop by adding the link ij and rejoining the loops, and [Ci, C 2 ] = otherwise. [Ci, C 2 ]ij 
denotes the loop that is made of C\ and C 2 with the addition of the link ij. When both C\ and 
C 2 are non-trivial loops, the third term describes a process where a loop splits into two loops (Fig. 
Q] (a)). When one of the two loops is a null loop, it describes a process where a loop becomes 
shorter by eliminating a self -retracting link (Fig. [U(b)). When both are null loops, it describes a 
self -retracting link disappearing (Fig. [0(c)). In the fourth term, Gij[C%, C 2 ] is a form factor that 
tells whether or not two loops C\ and C 2 are sharing the link ij : Gij [C±, C 2 ] = 1 if C\ and C 2 can 
be merged into one loop by removing the shared link ij, and G^fCi, C 2 ] = otherwise. (Ci,C 2 )ij 
denotes the loop that is made by merging C\ and C 2 by removing the shared link ij. The fourth 
term describes a process where two loops merge into one loop (Fig. [U(d)). In the small dz limit, 
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1/m 2 ~ 0{dz), and we can replace 4$ with in the third and fourth terms of the action to the 



'C 



linear order in dz. 

Note that double trace operators are generated for w^. Another set of auxiliary fields is intro- 
duced to express the double-trace operator in terms of single-trace operators as 



J c 



d(f) c du e 



-Ss 



(23) 



where 



S 5 = S uv [<f>P*,<f>P) 



N 2 



-N 



2m? 



Fij[Ci, C 2 ] 



\ci+C2]ij' 



(Cl+C2)ij 



+NM 2 tx{u\ j u ij ] 



c 



if) 



N 2 adz 
M 2 

-N^wc + NMhriul.Uij). 



(24) 



If we repeatedly apply the steps in Eqs. (fT3l) - (1241) to the last line of Eq. (|24l) i? times, we obtain 



n[< 

Z=0 



du e 



-S e 



(25) 



where 



adz 



c + aLcdztpQ 
C- 1 ) ^(0*^,(0* 



-A^g^ + jVM 2 ^^-)- (26) 

What is the physical meaning of the auxiliary fields ? In the last line of Eq. (|26l) . we note that 
0^ acts as a source for the low energy matrix field at scale e~ Rdz . The key difference from the 
standard renormalization group procedure is that the source fields are dynamical fields rather than 



fixed constants at each scale[12]. On the other hand, the equation of motion for <j>^' implies that 



< 



dR)* 
'c 



> 



N 



< Wc > ■ 



(27) 



12 



Therefore, the conjugate field <p ( c '* describes the Wilson loop operator. As we will see below, 
4>c and <f)* c are conjugate fields which satisfy a non-trivial commutation relation : sources and 
operators are conjugate to each other. 



(a) (b) (c) 

FIG. 2: Examples of the diagrams that contribute to the IR effective potential to the leading order in 1/N. 
Every link should be paired with another link with the opposite orientation. 



Finally, we integrate out u to obtain 

it 



n 

n[< 

1=0 



(0^(0*1 e -s 7> (2g) 



where 



S 7 = Suvi^AP] 
1=1 

+Si R [</> { < ? ) ]. (29) 
Here Sir is the effective potential given by 



Sir^} = -In J du e-^Mnviul^+N^c^ 



For a future use, we define 



(30) 



V[4$\ = ^s IR [^], (31) 



which can be computed using the strong coupling expansion, 

-l^UcUSM- ^ Lc >~5 Cl+ c 2+ c 3 ,o - .... (32) 
13 



Here the delta function is defined as 

= I I S QiAC],0i (33) 



<i,j> 



where Qij[C] is the U(l) charge defined on link ij associated with the flux of loop C [|13|l . If the 
loop C passes through the link ij from i to j (from j to i) n times, Qij[C] = n(—n). The first, 
second and third terms are from a self retracting loop (Figf2](a)), two loops (Figf2](b)) and three 
loops (Figf2](c)), respectively. Higher order terms can be obtained similarly. Now we take dz — > 
and R — > oo limits with (3 = Rdz fixed. Then, the partition function is written as 

D(p c D(j)* c e -( s buik&h( z )'0c(z)}+Suv{<f>h(o)^c(o)}+s IR i0c(P)]) ; (34) 
where 



Swuik = N 2 I dz 



'bulk 

a 



o 



M 2 



* c d z (j) c + aLc4>* c (j)c 
Fij 



(35) 



Since the partition function is independent of (3, we can take — > oo. From now on, we will 
interpret the scale parameter z as an imaginary 'time'. The dual description becomes a (D + 1)- 
dimensional field theory of closed loop. Although the action is written in terms of continuous z, 
one should go back to the discrete version whenever there is an ambiguity, e.g., when extracting 
boundary conditions by taking variations with respect to boundary fields. 

As is the case for matrix models, there are two important parameters that are independent with 
each other. The first is which controls the strength of quantum fluctuations of the loop fields : 
the whole action including the boundary actions scales as iV 2 . The second is the 't Hooft coupling. 
In this theory, there is no unique 't Hooft coupling. Instead there is a set of couplings defined 
in the space of loops, J{d,..,c n } which scales as the inverse of the 't Hooft coupling. Since we 
could have scaled out M 2 by redefining = U'^/M in Eq. ©, the theory depends only on 
the combination J{c 1 ,..,c n }^~'^" Lc ■ The small J{c lt ..,c„} limit is equivalent to the large M limit, 
which corresponds to the strong coupling limit of the matrix model where one expects to have the 
confinement phase. The set of J{c*i,..,c„}' s sets me magnitudes of loop fields in the bulk. We will 
see that background loop fields, in turn, control the size of strings which describe small fluctuations 
of the loop fields. 
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V. HAMILTONIAN PICTURE 



A. Partition function as a transition amplitude between many-body loop states 



The partition function can be viewed as an imaginary-time transition amplitude between many- 
body loop states. To see this, we will use a rescaled loop variable in this sub-section, 



$ c = n (pe- 



rn 



The bulk action in the new variable becomes 



''bulk 



-L 



dz 



$cd z $ c + aL c $* c <$> c 
NM2 



a 



(37) 



The action has the form for canonical bosonic fields, where $c ($c) corresponds to the coherent 
field associated with the annihilation (creation) operator defined in the space of closed loops. The 
annihilation and creation operators a c , cJ c satisfy the standard commutation relation 



ac, ay 



(38) 



where 5 C c > is a Kronecker-delta function defined in the space of loops. Then the partition function 
can be written as an imaginary-time transition amplitude, 



Z = lim < ^ffle'^l^i >, 
between the initial (UV) state at z = 0, 



(39) 



|*i> = j d<&* c d<& c V i [$* c ,<& c \\<& c >, 



(40) 



with 



(41) 



and the final (IR) state at z = oo, 



with 



I*/ > = J d§* c d§c y f l<5>* c ,®c}\®c >, 



(42) 



(43) 
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Here \l/i[$c> ®c] and tyf[&Q, $c] are the wavefunctions of loops written in the coherent state 
basis, 



|$ c >= e * ca c\0 >, (44) 

where |0 > is the vacuum in the Fock space of loops : a c \0 >= for all ac- (For the derivation 
of Eqs. (14TI) and (l43~T) . see Appendix. B). The bulk Hamiltonian is given by 



H = aL c ala c - j^^F ij [C 1 ,C 2 ]al i al; 2 a [ c 1+ c 2 ] t3 + G^Cu C 2 ]a\ Ci+C2) .. 0^0^.(45) 

The first term in the Hamiltonian describes a tension of closed loops. The second and the third 
terms are the interaction terms which describe the processes where one loop splits into two loops, 
and two loops merge into one loop, respectively, as is shown in Fig. Q] We use the convention 
a$ t = 0, aj = 1 for null loops. Similar loop Hamiltonians that describe joining and splitting 
processes of loops were considered in matrix models OOl 13 ill . 

This is an exact mapping between the D-dimensional matrix model {{D — 1) -dimensional 
quantum matrix model) and the (D + 1) -dimensional loop model (or £)-dimensional quantum 
loop model). Several remarks are in order. First, the Hamiltonian in Eq. (1431) is a many-body 
Hamiltonian that governs the quantum dynamics of loops along the scale z which is interpreted as 
an imaginary time. It is noted that the Hamiltonian is not Hermitian. Due to the cubic interaction 
term, the Hamiltonian is unbounded from below. However, the transition amplitude in Eq. (|39| ) 
is well defined because eigenvalues of the Hamiltonian are complex. Eigenvalues with a large 
negative real part in general come with a large imaginary part, and their contributions cancel 
with each other due to oscillation in phase. Second, the bulk Hamiltonian is universal, and it is 
independent of the details of the matrix model. All informations pertaining to the specifics of the 
matrix model are encoded in the initial wavefunction at z = 0. Third, the strength of the interaction 
between loops is order of 1/N, and loops are weakly interacting in the large N limit. Therefore, 
the theory becomes classical in the large N limit. Fourth, H does not have any hopping term such 
as a} Ci ac 2 with different C\ and CV This fact will become important for gauge symmetry, which 
will be discussed in Sec. V. For earlier works on string field theories formulated without quadratic 



action, see Ref. |32 



33] 



The fact that the partition function is independent of j3 has a remarkable consequence. By 
taking the derivative of Eq. (|39l) (for a finite (3) with respect to j3, we obtain 

= < 9 f \e- fiH H\^fi > . (46) 
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Since physical states are singlets of H, the Hamiltonian can be viewed as a generator of a 'gauge 
transformation'. The gauge transformation corresponds to a reparameterization of z. It is based on 
the fact that one can choose different speed of renormalization group flows at different scales with- 
out affecting the physics. By choosing the parameter a to be z-dependent, the reparameterization 
symmetry can be made explicit ni2fl . Here a(z) becomes the lapse function. Reparameterizations 
of z form a subgroup of the full diffeomorphism in the (D + 1) -dimensional space. It would be 
interesting to formulate the theory where the full diffeomorphism can be made explicit in the bulk. 
Here we proceed with the present formalism where we choose specific time slices along the z 
direction. 



B. Wilson loop operator 



sink 



z=0 



! 



FIG. 3: Loops are emitted at the UV boundary and propagate to the IR boundary. Since there is no hopping 
term in the Hamiltonian, loops can not move. Instead, they can join or split following the processes shown 
in Fig. [TJ Two loops with opposite orientations can get pair-annihilated through multiple interactions. 



The physical picture for the transition amplitude is the following. At the UV boundary (z = 0), 
a condensate of loops are emitted and propagate in z under the evolution governed by H. The 
amplitude of the condensate is < $ c < >~ O(N). This can be seen from the fact that the action for 
the unsealed loop fields <pc has N 2 as an overall prefactor, which implies < <pc >~ 0(1). Loops 
can join and split through the interactions as is illustrated in Fig. [3j A loop C and its anti-loop 
C, the loop with the opposite orientation, can get pair-annihilated through a series of interactions 
as is shown in Fig. |4](a). Moreover, a self -retracting loop can become a loop and an anti-loop as 
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(b) 



FIG. 4: (a) A process where a loop and its anti-loop get annihilated in pair through a series of interactions, 
(b) A self -retracting loop in the vacuum can split into a loop and an anti-loop. 

is shown in Fig. H](b). As it will be shown in Sec. VII. A, loop fields for self-retracting loops 
have non-zero vacuum expectation values in the bulk. Therefore, a pair of loop and anti-loop can 
be created out of vacuum. This means that two loops with the opposite orientations act as particle 
and anti-particle in a relativistic field theory. Finally, those loops emitted at the UV boundary are 
absorbed at the IR boundary. In this sense, the UV boundary is a source of loops, and the IR 
boundary is a sink. 



z=co 



z=0 




FIG. 5: A world sheet formed by multiple processes where a large loop absorbs many small loops at 
different stages to change its shape to disappear before it reach the IR boundary. 



Now let us consider a Wilson loop operator for a loop C which is much larger than the size of 
Wilson loops for which sources are turned on at the UV boundary. The expectation value of the 
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Wilson loop operator is given by the one-point function, 



1 d\nZ 

< Wr >= 

C N dJ c 



= lim — aiJ cl — . (47) 

j c =o < 9f\e-P H \% > 



If M is large, loops propagate independently in the bulk. To the zeroth order in 1 /M, the loop C 
propagate to the sink along the straight path. However, this configuration vanishes as e~ a ^ Lc in 
the large (3 limit because of the tension. In order for the expectation value to survive, the large 
loop C should absorb other smaller loops from the condensate to disappear before it reaches the 
IR boundary. Then the evolution of the Wilson loop forms a world- sheet in the bulk. One such 
configuration is shown in Fig. [5] Then the expectation value is given by the sum over all world- 
sheets of the Wilson loop. 

Since the interaction between loops is 0(1/N), loops become classical in the large N limit. 
This implies factorization of Wilson loop operators in the large iV limit, 

' n \ n 

W W c k ) =\{ (Wc k ) + 0(N n - 2 ). (48) 

\k=l I k=l 

VI. GAUGE SYMMETRY 

The absence of the hopping term in the Hamiltonian has a deep origin : the loop field theory has 
a gauge symmetry. Note that this gauge symmetry is not related to the U(N) gauge symmetry of 
the original matrix model. Loop fields are singlets for the U(N) gauge symmetry. In this section, 
we examine the consequences of the new gauge symmetry carefully. From now on, we return to 
the unsealed loop variable <j>c = ^ • 



The bulk action in Eq. (1351) is invariant under the time-independent transformation generated 
by Qij at each link 

<pc -»■ e ie » {i)Q " + » [c] <f) C , (49) 

where i is summed over all sites, /x is summed over D directions of nearest neighbor links, and 
9fj,(i) is a time-independent angle defined on the link < + >. The IR boundary action respects 
the symmetry, but the UV action does not. This is because the UV potential 



n=l {d,..,C n } 



n 



Ck 
k=l 



(50) 
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includes sources Jci,c 2 , -,c„ which explicitly break the symmetry. It is useful to view Jct,c 2 ,..,c n 
as an expectation value of another dynamical loop field. Then, the full theory is invariant if we 
allow the UV source to transform as 



J, 



Cl,C2,--,C n 



e ^(i)ELiQtf+4 c ).] j 



Ci,C2,--,C n - 



(51) 



This time-independent symmetry can be lifted to a full space-time gauge symmetry by introduc- 
ing temporal components of a two-form gauge field B MN in the bulk with M, N — z, 1, 2, D, 



?bulk 



N 2 j dz (d z + iQ ll+ ^C}B^ z (i, z) S j(t> c + aL 



c<P c <Pc 



a 



M 2 



Fij [Ci,C 2 ]0 

c 1 4 > c2 ( t ) [c 1 +C2\ l j + Gij[Ci, C 2 ](t)\ci+C2)ij ■ ( t ) Ci ( t ) C2 



(52) 



where B^ z {i) with /i = 1, 2, .., D are the temporal components of the two-form gauge field defined 
at each spatial link. This two-form gauge field is the Kalb-Ramond gauge field ll29n . Now the full 
theory is invariant under the space-time dependent gauge transformation with 

Bfufaz) ->• B^ z (i,z)+ 9 z (i + fi,z)-6 t (i,z) - (d^i, z)), (53) 

where 9 z (i, z) is a temporal gauge parameter defined at each site. This is the discrete version of 
the usual gauge transformation for the two-form field, Bmn Bmn + 9m9n — 9n0m- Note 
that introducing the temporal components of the two-form gauge field into the theory doesn't do 
anything except for making the gauge symmetry more explicit. This can be understood from the 
fact that one can reproduce the original action in Eq. (|35l) by choosing the temporal gauge with 
Bfj. z — 0- This can be done by choosing 



dz Bfuii^z 



(54) 



with 9 z (i, z) = in Eq. (|53l) . The temporal components can be completely gauged away because 
they are pure gauge degrees of freedom in the presence of boundaries. This is in contrast to the 
case with the periodic boundary condition, where the time independent component of the temporal 
gauge field can not be gauged away. 

As a result of the gauge symmetry, there is no quadratic hopping term for loops in the Hamil- 
tonian. However, this does not necessarily mean that loops are always localized in space. Loops 
can change their shapes and move in space by absorbing or emitting other loops. For example, 



20 



FIG. 6: Loops can change their shapes and move in space like amoebas, by absorbing or emitting small 
loops. 



Fig. [6] shows a loop changing its shape by absorbing two small loops. Therefore, loops can prop- 
agate with the help of other loops. If loop fields are 'condensed', whose precise meaning will 
become clear in a moment, the condensate provides a coherent background on which other loops 
can propagate. Loops propagate 'on the shoulders of other loops' to explore the bulk space. This 
is analogous to the the slave-particle theory discussed in Sec. II. One difference is that loop fields 
themselves play the role of 'hopping fields' for other loops, while in slave-particle theory the hop- 
ping field is a bi-linear of slave-particle fields. The difference originates from the fact that loops 
are extended objects while slave-particles are point objects. Only when the condensates of loop 
fields are 'coherent', the bulk space is regarded as a well defined extended space by loops. Oth- 
erwise, loops are more or less localized in space. In this sense, an extended space emerges in the 
bulk as a dynamical feature of a phase where loop fields form coherent condensates. 

When do loop fields become coherent ? To make this notion more precise, we first note that 
the phase modes of complex loop fields <pc = \4>c\e lbc play the role of the spatial components of 
the two-form gauge field. To see this, suppose that the loop field has a background value < <p c >• 
Then the cubic interaction generates a quadratic hopping term, 

M^ a c+c' a C- ( 55 ) 

The amplitude of < <p c > is the strength of the hopping, and the phase determines the geometric 
phase acquired when the loop C' hops to C + C. Therefore be plays the role of the spatial 
components of the two-form gauge field to which loops are electrically coupled. Note that the 
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two-form gauge field is also a part of dynamical loop fields. We identify 

bc= f B, (56) 

J A c 

where B^ u is the spatial components of the two-form gauge field and the integration is over an 
area Ac enclosed by the loop C . Let us focus on the loops with unit plaquettes in which case we 
take Ac as the surface spanned by the unit plaquette. 

Although the two-form gauge field does not have the bare action, it acquires the kinetic energy 
from quantum fluctuations. This is similar to the way that the Maxwell's term is dynamically 
generated for the auxiliary gauge field in the slave-particle theory as discussed in Sec. II. The 
gauge coupling for the two-form gauge field is renormalized to 0(1/N 2 ). This can be understood 
by integrating out 'heavy' loop fields to obtain an effective action for 'light' loop fields in the bulk. 
It is easiest to see the generation of the kinetic energy in the large M limit, where we can use 1/M 
as an expansion parameter. The 'mass' of a loop field is proportional to the length of the loop 
because of the tension. We integrate out loop fields with L > 4 and obtain an action for the loop 
fields with L < 4. In particular, we focus on the effective action for the shortest non- self -retracting 
loops with L c = 4 whose phase modes can be viewed as the spatial components of the two-form 
gauge field on unit plaquettes. For simplicity, we choose the temporal gauge with B^ z = and the 
scale of z to set a = 1. 

Let us consider three vertices in Eq. (|52l) . Each vertex has the form ^§c i §c il §Ci2 with 
% — 1, 2, 3, where Cj's have length 6 and Cu's with I — 1, 2 have length 4 as is represented in 
Fig. [7] (a). They describe the processes where a loop on a unit plaquette with sides //, v merges 
with a loop on a unit plaquette with sides v, A to form a loop with with length 6. Here we interpret 
0Ci' s as heavy fields and 0q ; ' s with I — 1, 2 as light fields. In particular, the phase modes of <f>c a 
represents the two-form gauge field defined on each plaquette. Now we integrate out the heavy 
loop fields using the quadratic action. Because this quadratic action has the local U(l) symmetry 
in the loop space, 0c — > e^cpc, we need to introduce a series of vertices in order to saturate <p c 
with <pc and obtain a non- vanishing result. A minimum path to saturate all heavy fields is shown in 
Fig. 13(b). In the first step, we add a vertex of the type |^0* Ci+C2 ) 4j 0Ci0c 2 > where {C\ + C 2 )y is 
a loop that results from merging G\ and C 2 by removing one shared bond. Integrating out (pc 1 and 
0c* 2 , we obtain the loop fields in the second configuration in Fig. |7](b). In the second step, we use 
a vertex ^^* 0l+ c 2 )'. 0(<7n-cb)«> and integrate out (j)( Cl +c 2 ) ir In this step, the merged loop in the 
first step become a shorter loop (C\ + C 2 )'ij by eliminating one self-retracting link. The remaining 
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(a) 




(b) 

FIG. 7: Kinetic energy for the two-form gauge field generated by heavy loop fields, (a) Three vertices each 
of which involves one loop with length 6 and two loops with length 4 can generate a term for the six loops 
with length 4, once the loop fields with length 6 are integrated out. The resulting term for the small loops 
becomes the kinetic energy for the two-form gauge field, (b) In integrating out the loop fields with length 6 
in (a), one has to introduce a series of nine vertices. Each step depicts a process of adding a new vertex and 
integrating out one or two loop fields : two fields for the first and the fourth steps, and one loop field for the 
other steps. 

steps can be understood in a similar way. In total, nine vertices and eleven loop propagators are 
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needed. (Note that each of the first and fourth steps introduces two propagators because two loop 
fields are integrated out in those steps, while all the others involve only one propagator.) Each 
vertex contributes |p and each propagator contributes Combined with the factor (^j^j from 
the original three vertices, we obtain an action for the light loop fields, 



S e ff ~ — / dz 



N 2 



M 24 
cubes 



i=l i=l 



(57) 



where the summation is over all cubes in the D-dimensional lattice. The second term is from the 
same process for the anti-loops. The loop field for the anti-loop, 0c is in priori independent of 0c- 
However they are dynamically mixed. Because of pair-annihilation and pair-creation processes of 
loops and anti-loops as is shown in Fig. HI the effective action should include terms of the form, 
0c0c an d fic^h- ^ s a result' me phase modes of 0c and 0c are locked. Only the anti- symmetric 
mode with b c = —be remains gapless in the presence of mixing. If 0c«'s have finite amplitude 
0o, this gives the standard 'magnetic' term for the two-form gauge field defined on each cube of 
the lattice 

_ J_ f dz cos [« 3 ( A m^a + A„B Xfl + A X B^)] , (58) 
9kr 3 cubes 

where we use the fact that b D = a^B^ for a unit plaquette with sides ji, v. The finite derivative is 
defined as = b vX {x+x*)-b vX {x) ^ Rere ^ ^ j^4>t is the renormalized coupling for 

the Kalb-Ramond (KR) two-form gauge field. 

Now we turn our attention to the 'electric' term which involves the time-derivative of the gauge 
field. We consider the quadratic action for the loop fields on unit plaquettes, 

N 2 [4>h(d z + L c )4>c + 4>h(d* + L c)<Pc] ■ (59) 

If we integrate out the amplitude fluctuations of the loop fields, the time derivative term will be 
generated for be- Because of the dynamical constraint b c = —be caused by the mixing between 
0c and 0c, the linear time derivative term is canceled, leading to the second derivative term, 

iVV I dzJ2(d z B^) 2 (60) 

for each plaquette. Eqs. (l58l) and (l60l) represent the full kinetic energy term for the two-form 
gauge field in the temporal gauge. 

Because of the gauge symmetry, the mass term is not allowed for the two-form gauge field 
in the bulk. However, the gauge symmetry does not automatically imply that the two-form field 
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(D+l)-dimensional bulk 



NS D4 -brane 



FIG. 8: An A r 5o-4-brane in the (D+l)-dimensional space. Around an A r S*£)_4-brane, there is a net flux of 
2-7T for the three-form flux : J s3 H = 2it. When the tension of the AS^^-brane is positive, A5D_4-branes 
generated out of quantum fluctuations remain small. When the tension becomes negative, A£D_4-branes 
are proliferated in the bulk. 

arises as a massless excitation in the bulk. This is because of the compactness of the phase mode 
: b c ~ be + 2ir, which allows for a topological defect to exist as a magnetic excitation of the 
gauge field. In the presence of topological defects, the field strength H = dB does not satisfy 
the Bianchi identity dH = 0. In (D + 1) -dimensional space-time, the topological defect which 
carries a magnetic charge is a (D — 4)-brane, which is a (D — 3) -dimensional object in space-time. 
We call this object NSn-4 brane ll52ll . Around the A r S , £ ) -4-brane, there is a net 2n flux for the 
three-form flux, 



dH(x) = 2tt / d£ D - 3 5 D+1 (x - 0, 



(61) 



where £ is the coordinate of the NS-brane embedded in the (D + 1) -dimensional space, and d£ D ~ 3 
is the oriented volume element of the brane. This is illustrated in Fig. [8] Note that the Dirac 
quantization condition between the charge carried by loop fields, which is set to be 1 as can be 
seen from Eq. (|49l) , and the charge of the A^S^^-brane is automatically satisfied. This follows 
from the fact that the phase 2n on a unit plaquette is invisible to loop fields. In D > 4, this is a 
brane extended along (D — 4)-directions in space at a given time slice with fixed z. In D = 4, this 
is a point-like particle. In D = 3, this is an instanton which is localized both in space and time. 
Whether loop fields provide a coherent background for other loops is determined by dynamics of 
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A^S , D-4-branes. 




(a)z<z NS (b)z = z NS (c)z>z NS 



FIG. 9: An iVSD-4-brane can be viewed as an instanton in D = 3 where one inserts a source of 2ir 
magnetic charge to the theory at a given z^s- ( a ) For z < z^s, the phases of loop fields are zero, (b) 
At z = ztvs> the phases of loop fields are 2ir along a semi-infinite line creating a cube which contains the 
3-form flux of 2tt : H123 = A1.B23 + A2-B31 + A3B12 = ^S x ,x NS - ( c ) F° r z > z^s the configuration in 
(b) is smoothly deformed and the flux is smeared over a region which contains the net flux 2-k. The size of 
red dots represents the amount of 3-form flux contained in each cube. 

The physical nature of A^S , £ ) _4-brane can be most easily understood in D = 3 where NSo-i- 
brane is an instanton localized at a point in the four-dimensional bulk. We start with a configuration 
of the loop fields D fof unit plaquettes with 

b ni2 (x, z) = 2tt5 x i )X i ns 5 x 2 :X 2 ns Q(x 3 - x 3 NS )Q(z - z NS ), (62) 

where the phases of the loop fields on x 1 — x 2 plaquettes are 2it along a semi-infinite line in the 
3-dimensional space for z > z^s, and the phases are zero, otherwise. Since be = be + 2n, this 
configuration is equivalent to the trivial configuration where be = everywhere. However, we 
can view this configuration as a topological defect with a non-trivial three-form flux on a cube, 
#123 = 2tc5 3 (x — xns) © (z — zns)- This means that there is a source of magnetic flux for the 
two-form gauge field localized at a point in the bulk, dH = 2it5 3 (x — x NS )8(z — Zns)- Since 2ir- 
magnetic flux is concentrated at one cube, the topological defect is trivial. Now the configuration 
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n = 2% 

s 3 



z=0 





FIG. 10: A pair of instanton and anti-instanton in the four-dimensional bulk for D = 3. 

is deformed smoothly, b' c = be + Sbc- Under a smooth deformation, the flux is smeared out over 
an extended region in the space, while the net flux 2% does not change. Now the flux is visible 
by loop fields. This is illustrated in Fig. [9] As z is increased further, the flux can merge back 
into one cube and disappear into the vacuum through the inverse process. This describes a pair 
of instanton and anti-instanton as is shown in Fig. [TOl In higher dimensions, A^S'D-4-branes are 
extended objects. In (D + 1) -dimensional bulk, we can think of Fig. [10] as a configuration in a 
slice at a fixed x 5 , ...,x D+1 , where the NSn^-brane is extended along the (D — 3) directions. 
They can be wrapped into compact objects as is shown in Fig. [8l which in a sense describe bound 
states of A^S'i5_4-brane and anti-A^S , D_4-brane. If the size of the wrapped A^S'D-4-branes become 
infinite, A^S^^-brane and anti-A^5 , D-4branes become unbound. 

The tension of the A^S^^-brane is proportional to A^ 2 because A^S , D-4-brane is a topological 
defect of the two-form gauge field with the coupling proportional to 1/N 2 . Here we are using the 
term 'tension' in a loose sense. In D > 4, it literally means the tension of A^S , D_ 4 -branes. In 
D = 4, it refers to the mass of W5*£i_4-particle'. In D = 3, it refers to the action of l NSn-4- 
instanton'. For a sufficiently large N, we expect that A r S , £>_4-branes are gapped. In this case, 
A r S'i5_4-branes will be wrapped into compact objects with a finite size in the vacuum. For a 
small N, the bare tension of NSn-i-brane is small, and quantum fluctuations can renormalize the 
tension into a negative value. Then A^S , £ ) _4-branes are condensed, and extended A^S^^-branes 
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fill the space in the bulk. It is also possible that iVS£)_4-branes always condense for any finite N 
in low dimensions. We will discuss the consequences of different dynamics of A^5*D-4-branes in 
the following section. 



VII. EMERGENT SPACE AND QUANTUM ORDER IN HOLOGRAPHIC PHASES 

In this section, we will discuss various phases that the matrix model can have, by focusing on 
the behavior of iVS^zt-branes. In particular, we will see that the dynamics of string excitations 
around a saddle-point configuration of the loop fields is determined by the fate of A^S , £)_ 4 -branes 
in the bulk. In order to discuss about this issue systematically, we first turn to the saddle point 
equations. 



A. Saddle point solution 



The saddle point configuration of loop fields is determined from the equation of motion. 



d^c = -2, C C + — ^-[^C^^ (63) 
-d z <j> c = -Wc+jlp {FiAOuC^^ .(64) 

These equations are supplemented by two sets of boundary conditions. It is more convenient to use 
the action with discrete time step dz to isolate boundary fields from bulk fields. The UV boundary 
condition is obtained from Eq. (fT2l) . 



dSx 
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(65) 



and the IR boundary condition from Eq. (|29l) , 

dS 7 



d4 R) 



N 



J c 



c W } 



0. 



(66) 



When the UV potential includes only single-trace operators, V = —Jc<fic , Eq- (1651) leads to the 
standard Dirichlet boundary condition for the source field : 4>q = Jc- For more general non-linear 
UV potential, it becomes a mixed boundary condition. This is consistent with the prescription for 
the UV boundary condition in the presence of multi-trace deformations in the standard AdS/CFT 
dictionary [ | 34f | . The IR boundary condition is a mixed one because V' is in general non-linear. 
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For self-retracting loops, V also contains terms that are linear in loop fields as is shown in the 
first term in Eq. (l32l) . Eq. (|66|) then implies that <p* c ^ at the IR boundary for self-retracting 
loops. As will be shown in the next paragraph, this means that loop fields for self -retracting loops 
have non-zero expectation values at all z in the bulk. This, in turn, generates non-zero vacuum 
expectation values of the source fields <\> c for self-retracting loops. As was discussed in Fig. HI 
self-retracting loops can turn into a loop/anti-loop pair through an interaction. 

In general, the saddle point configuration is ^-dependent, and <\> c (z) is not necessarily the com- 
plex conjugate of 4>* c (z) . One should treat 0c and C as two independent fields. Then the equations 
of motion can be viewed as a set of Hamiltonian equations in the phase space of {0c, 4>c}- 

A 

- ~J 



~J/M' 







FIG. 11: A schematic profile of a loop field 0c and the conjugate field C in a deep confinement phase. 
0c and 4? c satisfy the boundary condition given by Eq. (l67l) at any z in the bulk. 



Although Eqs. (1631) and (|64l) are coupled equations for c and C , one can eliminate <p* c in favor 
of 0c- We first note that the partition function and all observables including vacuum expectation 
values of Wilson loop operators represented by 0c(O) are independent of how we choose (3 in Eqs. 
(|34|) and (1351) . This means that the saddle point solution (j>c{z) and 4>* c (z) for z < (3 is independent 
of /3. Since we could have put /3 anywhere, 0c* and <p* c should satisfy the IR boundary condition 
at any z, 

dV'[<t> c {z)\ 



Fob) 



» 



(67) 



This is illustrated in Fig. [HJ The fact that one can put the IR boundary at any z has an interesting 
implication on the role of the IR boundary. Usually, one can associate a boundary condition with a 
physical object located at the boundary. However, Eq. ([67]) is special in the sense that an observer 
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at z < (3 can not 'feel' the presence of a physically identifiable object at z = (3. Suppose one stops 
the renormalization group procedure at z — f3 and impose Eq. (1671) at the IR boundary. If a UV 
observer sends a wave toward the IR region, the reflected wave from the IR region is exactly the 
same as the reflected wave one would observe in the space which is extended to z = oo without 
an boundary. In this sense, the IR boundary is not a physical boundary : one can always trade the 
IR boundary with the space where z is extended to infinity. 

Using Eq. (1671) . we can write a set of first order differential equations for the source field only, 



d z <p c = -L c <Pc + ^ (-2F tJ [C i c /^ Cl V [c + c l]iJ + G^C^c^cAc^c^ ■ 

(68) 

Once <pc (z) is solved using the UV boundary condition in Eq. (1651) . the conjugate field is readily 
determined from Eq. (1671) . One can check that the source and the conjugate field satisfy Eq. (1671) 
at all z through an explicit calculation perturbatively in 1/M (see Appendix C). 

It is tempting to interpret Eq. (|68l) as the beta function of the sources for Wilson loop operators. 
However, there is an important caveat for this interpretation, which comes from the fact that this 
is the saddle point equation of the quantum theory for dynamical loop fields. The saddle point 
equation is expected to be valid only when quantum fluctuations are weak for a sufficiently large 



N. For a smal 



N, one can still have a well-defined beta function under the usual renormalization 



group flow 11271 12811 . However, the beta function can not be directly identified with the saddle 
point equation of the loop fields if the saddle point solution becomes unstable by strong quantum 
fluctuations. As we will see in the following sections, non-perturbative fluctuations can invalidate 
the holographic description for small N. 

B. Fluctuations near the saddle point 

Fluctuations near the saddle point configuration <p c (z) describes dynamical string in the bulk, 

(f>c(z) = te(z)+Xc(z), (69) 
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where xc describes small fluctuations around the saddle point. We call xc string field to distin- 
guish it from the loop field 4> c . The dynamics of string is governed by the action, 



Sbuik = N 2 [ dz x* c dzXc + LcXcXc 
Jo l 



■— (F i:i [Ci, C 2 ](p [ c 1 +C2h J ( z )x* Cl X*c 2 + 2F ij [C u C 2 ]4>c 1 ( z )Xlc 1 +c 2 ] i Xc 2 
+ Gij[C x , C 2 }4>\ Ci+C2hj (z)x Cl Xc 2 + 2Gy[Ci, C 2 ]0 Cl (z)xc 2 X*c 1+ c 2 x 

1 



FijlCx, C 2 \x[c 1 +c 2 ] ij Xc 1 Xc 2 + Gij[C u C 2 \x Cl Xc 2 x\c^c 2 ) % 



(70) 



M 2 

Here Xc{z) = Xc( z ) + ^Xc( z ) ls a complex field. Following the standard method of the steep- 
est descent, the contours of the real and imaginary parts of the complex fields are chosen so that 
the real part of the eigenvalues of the quadratic action becomes maximum along the deformed 
contours p^i . Note that the string fields xc acquire the hopping term through non-zero con- 
densates of loop fields. It also has the terms that describe pair creation/annihilation of two closed 
strings. 



C. Possible phases of the matrix model 

In this section, we describe possible states of the matrix model using the holographic descrip- 
tion. In particular, we will see that strings have different dynamics depending on the behavior of 
iVS , D_ 4 -branes. One observable that is useful in distinguishing different states is the correlation 
function between Wilson loop operators. In particular, we focus on the correlation function of 
phase fluctuations of Wilson loop operators, 

F(C,C') = (5b c 5b c ,), (71) 

where 5bc = be— < be >■ In the bulk description, this correlation function corresponds to a 
two-point string-string correlation function. This object is of particular interest because the string 
state that corresponds to the phase mode describes the two-form gauge field in the bulk. 



1. Confinement phase 

For non- self -retracting loops, V' is quadratic or of higher order in the loop fields. For small 
</>c(0), the first term on the right hand side of Eq. (1681) dominates. As a result, <j>c(z) decays 
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(c) (d) 



FIG. 12: The Wilson loop correlation function in different phases, (a) In the confinement phase, NSd-a- 
branes are condensed in the bulk, and the string emitted at the boundary stays near the boundary. Due to 
weak fluctuations of string world sheet, the string propagates in a straight path, (b) In the non-holographic 
critical phase, iV5£>-4-branes remain condensed in the IR region. However, large amplitudes of loop fields 
in the UV region make the string world sheet highly fluctuating along the D dimensions. As a result, the 
correlation function decays algebraically, (c) In the holographic critical phase, loop fields acquire finite 
expectation values everywhere in the bulk, and A r S , £j_4-brane is gapped. Strings can propagate deep inside 
the bulk, and the correlation function shows a power-law behavior with a scaling dimension determined 
by the mass of the string, (d) In the deconfinement phase, loop fields with infinitely large loops acquire 
non-zero expectation value in the bulk. This causes highly non-local fluctuations of the world sheet of a 
string in the bulk. 

exponentially in z. Because larger loops for which sources are not turned on at the UV boundary 
are generated out of many small loops, amplitudes with larger loops decay exponentially with the 



32 



area enclosed by the loop. This corresponds to the confinement phase of the matrix model. In the 
confinement phase, one can define a cross-over scale z* beyond which loop fields have negligible 
amplitudes. 

Now let us consider dynamics of strings in the IR (z > z*) and the UV (z < z*) regions 
separately. In the deep IR region, amplitudes of non- self -retracting loop fields are exponentially 
small. This means that loops that are emitted from the UV boundary rarely reach the IR region. In 
this region, phase fluctuations of loop fields do not have a stiffness, which leads to a condensation 
of iVS , D_4-branes. As a result, strings are subject to the strongly fluctuating two-form gauge 
field. In this region, strings are confined, and the two-form gauge field is gapped l!37l 13811 . For 
a more systematic discussion on possible phases of anti-symmetric fields, see Ref. [39]. In the 
deep IR region, strings can not propagate by themselves; only charge neutral bound states of string 
and anti-string can propagate. On the other hand, loop fields have significant amplitudes in the UV 
region. The source field J{c lr .,c„} plays the role of a symmetry-breaking field at the UV boundary. 
As a result, phase fluctuations of loop fields are small, and iVS^/i-branes are suppressed in the 
UV region. Because loop fields are coherent near the UV boundary, strings are deconfined in this 
region. There is a domain wall that separates the IR region with condensed iVS , D_4-brane and the 
UV region without iV£D_ 4 -brane. 

In the confinement phase, a string that is emitted from the boundary can not penetrate through 
the wall of condensed iVS/j^-brane. Therefore it stays within the UV region. Deep inside the 
confinement phase with small J{Ci,..,c n }, the condensates of loop fields are small. As a result, the 
hopping amplitudes of strings are small, and fluctuations of string world sheet is small. For the 
correlation function in Eq. (fTTT) . the strings inserted at the UV boundary are connected through a 
minimum number of hoppings, forming a straight path as is shown in Fig. [[21(a). This leads to 
an exponentially decaying correlation function for the Wilson loop operators. In the confinement 
phase, the bulk geometry ends at a finite scale due to the proliferation of iVS^^-branes. This is 
reminiscent of the idea that geometry can get truncated by tachyon condensation[40]. 

As J{d,..,c n } is dialed up, amplitudes of loop fields in the bulk increase. Accordingly the 
cross-over scale z* increases. At the same time, fluctuations of string world sheet increase as 
the amplitudes of loop fields become larger. Suppose the system becomes critical either by fine 
tuning or dynamical tuning. In the case of fine tuning, one may have to tune more than one 
microscopic parameters to reach a critical point. For the following discussions which focus on 
physical properties of the critical states, it is not important whether those states are realized as 
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phases or critical points. So we will use the term 'critical phase' in a broad sense to include 
not only critical phases realized within a finite region in the parameter space of a microscopic 
model but also critical states realized at critical points by fine tuning. Logically, there exist at 
least two different scenarios via a criticality is achieved. In the first scenario, iVS^-^-branes 
remain condensed in the IR region with a finite z*. But loop fields acquire large amplitudes in the 
UV region so that strings are delocalized along the D-directions, mediating critical correlations 
between operators inserted on the UV boundary. In the second scenario, the cross-over scale 
diverges and A^S , D_ 4 -branes are suppressed out in the bulk. In this case, strings can propagate 
deep inside the bulk. In the following, we will discuss the two scenarios in more detail. 



2. Non-holographic critical phase 

Here we discuss the first case which is likely to be realized when N is small. For a small N, 
iVS£>_4-branes are 'light'. Even though loop fields have finite amplitudes in the bulk, quantum 
fluctuations may destabilize the saddle-point solution. Indeed this is what always occurs in the 
pure 2-form gauge theory in the flat four dimensional space|37, 41, 42]. If the mass or tension of 
iVS , D_4-brane becomes negative, iVSz^-branes are condensed. Once iVS/^-branes are prolif- 
erated, the two-form gauge field acquires a mass gap, and strings are confined in the bulk, as is the 
case in the confinement phase. The difference from the confinement phase is that the boundary 
theory is critical. The boundary theory can be critical although strings are confined deep inside the 
bulk because critical correlation between boundary fields can be mediated by strings that propagate 
near the UV region. It is noted that 'confinement' in the boundary matrix model and 'confinement' 
of strings in the bulk are not the same thing. In this phase, strings are confined in the bulk, but the 
matrix model is not in the confinement phase. We call this non-holographic critical phase. 

In this critical phase, a string emitted from the UV boundary no longer takes the straight path 
because of large amplitudes of background loop fields in the UV region. Rather, the world sheet 
of string strongly fluctuates, and the correlation function can decay in a power law because of 
delocalized strings. However, strings are still localized within the UV region along the z direction 
due to the condensation of iVSD-4-branes in the IR region. The fluctuations of the world sheets of 
strings is predominantly along the D space dimensions as is illustrated in Fig. fT2l(b). 
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3. Holographic critical phase I : deconfined closed string 



The second scenario is qualitatively different from the previous one. In this phase, loop fields 
develop non-zero amplitudes in the IR region. For a sufficiently large N, iVS£>_4-branes are 
suppressed, and the saddle-point solution is stable against non-perturbative fluctuations of the 
two-form gauge field. If iVS£>_4-branes are suppressed, the compactness of the gauge field is 
unimportant at long distances. Strings are deconfined in the bulk because loop fields provide a 
background in which strings can propagate coherently. Note that strings are still coupled with the 
dynamical two-form gauge field, but the gauge field is no longer confining in this phase. This is 
analogous to the Coulomb phase discussed in Sec. II. Thanks to the coherent background loop 
fields, closed strings can explore the extended space in the bulk. The bulk space is not a gauge 
invariant object. However, it assumes a 'classical identity' in the large N limit where fluctuations 
of loop fields are suppressed. In this sense, the bulk space emerges in the holographic phase, but 
not in the non-holographic phase. 

Because of the source that explicitly breaks the gauge symmetry, the transformation B MN — > 
Bmn + 9m9n ~ 9^9m, is not a symmetry at the UV boundary. As a result, the U(l) mode 
9 M becomes a physical mode. This means that there is a U(l) gauge field localized at the UV 
boundary in the dual theory. This £7(1) mode originates from the Abelian component of the U (N) 



44, 



matrix theory. We identify this as the singleton mode localized at the boundary [7L 

It is noted that the symmetry breaking source J{c 1 ....c n } at me UV boundary does not neces- 
sarily open up a gap for the phase modes of loop fields in the bulk. This is because the symmetry 
breaking source is only at the UV boundary, but not in the bulk. This is analogous to the case 
where one applies a symmetry breaking field at the boundary of a system where a global sym- 
metry is spontaneously broken in the bulk. Although the boundary field determines the direction 
of the symmetry breaking in the whole system, the Goldstone mode in the bulk survives in the 
thermodynamic limit. 

The dynamics of the D-dimensional matrix model in the long distance limit is governed by 
strings that propagate in the (D + 1) -dimensional space. We call this phase holographic phase. 
The hallmarks of the holographic phase are deconfined strings that propagate in the bulk with 
the extra dimension and the emergence of the Bianchi identity dH = for the two-form gauge 
field in the long distance limit. It is emphasized that these features are not protected by any 
symmetry of the original matrix model. They are dynamical properties which emerge only in 
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the holographic phase. The quantum order present in the holographic phase is analogous to the 
quantum order associated with the emergent Bianchi identity in the fractionalized phase of the 
slave-particle theory as discussed in Sec. n. 

In the holographic critical phase, strings can propagate deep inside the bulk as is shown in 
Fig. [121(c). The correlation function shows a power-law decay through a classical trajectory that 
is extended to the bulk. Because of the gauge symmetry and the non-trivial quantum order, we 
expect that the scaling dimension of the phase mode of Wilson loop operators will be protected 
accordingly. To determine the scaling dimension, one has to first solve the loop equations in the 
bulk and find the string Green's function in the background determined by the loop fields. We 
defer an explicit calculation for future studies. 

In the four-dimensional J\f = 4 super Yang-Mills theory, the phase mode of Wilson loop oper- 
ators has the scaling dimension 6 ll47n . The reason why it is not 2, which is the expected scaling 
dimension for massless two-form fields in D = 4||7|], is that a Chern-Simons coupling generates a 
mass for the two-form gauge field through the mixing with the Ramond-Ramond fields. However, 
the scaling dimension is still protected from acquiring a large quantum correction in the large N 
limit. Such protection of scaling dimension is often attributed to supersymmetry. However, the 
non-trivial quantum order will protect the scaling dimension of the phase mode from acquiring a 
large quantum correction even in non- super symmetric holographic critical phases in the large N 
limit. This is true whether or not the two-form gauge field becomes massive through a Chern- 
Simons coupling with Ramond-Ramond fields. This is because the string theory becomes classi- 
cal in the large N limit, and the coefficient of the Chern-Simons term is quantized. Therefore the 
mass of the two-form gauge field can not become large even when other string modes become very 
massive at strong coupling of the boundary matrix model. This, in turn, implies that the scaling 
dimension remains small for the phase mode of Wilson loop operators. This is in sharp contrast to 
the non-holographic critical phase where it is expected that the operator generally receives a large 
quantum correction at strong coupling. 



4. Deconfinement phase 

Strictly speaking, the critical phases discussed in the previous two sections are kinds of decon- 
finement phases. Here we use the term 'deconfinement phase' in a narrower meaning, that is, free 
theory in the IR limit. If the sources at the boundary are very large, the second and third terms 
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in Eq. (1681) dominate, and the source fields will grow as z increases for D > 4 in which case the 
boundary matrix model is expected to flow into IR free gauge theory in the weak coupling (large 
J{Ci,..,c n }) limit. As the amplitudes of the loop fields become larger, large loops are generated 
through the joining processes. As a result, loop fields with all sizes are condensed in the bulk. 
Then strings propagating in the bulk become highly non-local because strings can hop from one 
configuration into another configuration which is very different from the initial one. This is a 
string condensed phase. In this phase, the two-form gauge field acquires a mass and iVSo-^brane 



is confined due to the Higgs mechanism [|43L 



46]. 



In the deconfinement phase, a string emitted from the boundary becomes very large in the bulk 
and lose its identity as a closed string. The critical fluctuations are mediated by highly non-local 
fluctuations in the bulk. In this phase, the locality is lost in the bulk. 

The deconfinement phase can be viewed as an extreme limit of the holographic critical phase 
discussed in the previous section. Even in the holographic critical phase, some loop fields are 
condensed in the bulk, as is the case in the deconfinement phase. The difference is that only small 
loops are condensed in the holographic critical phase while loops with all sizes are condensed in 
the deconfinement phase. Note that condensations of small loops do not generate a mass gap for 
the two-form gauge field. This is because closed strings with finite sizes as point-like particles 
are coupled only with the field strength tensor of the two-form gauge field. In certain models, 
one can in principle change microscopic parameters to tune the size of condensed loops, smoothly 
interpolating between the holographic critical phase and the deconfinement phase. This basically 
controls the size of strings in the bulk. The N = 4 SU(N) gauge theory in four dimensions is 
believed to be in this class : for a sufficiently large N, one can smoothly tune the 't Hooft coupling 
A from a large value to zero without going through a phase transition. The one parameter family 
of the critical theories form a line of fixed points. Here, A = is a special point where the size 
of string diverges. In non-supersymmetric theories, it is expected to be harder to stabilize a theory 
at an arbitrary gauge coupling. Most likely, we expect that the holographic critical phase will 
arise as a multi-critical point between the confinement phase and the deconfinement phase for a 
sufficiently large N. 
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D. A mean-field description 



Some features of the phases discussed in the previous section can be easily understood if we 



focus on a subspace within the space of loop fields {4>c}- We focus on the mean-field Ansatz[43, 



460 where a loop field is represented by a product of link fields along the loop, 



<Pc{z)= J] ^(i,z). (72) 

(i,i+n)ec 

Here £ M (z, z) = + z) is a complex scalar field defined on the link i,i + //. This is a 

huge simplification where we reduce the space spanned by functions defined on loop space into 
the space spanned by functions defined on the links. Under the gauge transformation, the link 
variables transform as 

£ M (i,z)^e*^>e M (i,s). (73) 

Therefore, these link fields should be charged with respect to the two-form gauge field. A minimal 
action for the link field that has the same symmetry as the original loop model is an Abelian-Higgs 



model[43] for the link field, 



S b uik = -tJ2^M(i)^N(i + M)e M (i + N)e N (i)e- iBM ^+c.c.] + £V(|&if(i)| 2 ) 



i,M 



— ^2 cos(A L B MN + A M B NL + A N B LM ). 



(74) 

9kr cubes 

Here we discretize the z direction, and the action is written in a (D + 1) -dimensional lattice. The 
box represents sum over all plaquettes including temporal plaquettes. The link fields £ z along the 
temporal directions can be viewed as an auxiliary field that is introduced to keep the two-form 



gauge symmetry. One can use a continuum description as well[46]. 

The phase structure of this model is very similar to the one for the Abelian-Higgs model for 
scalar fields discussed in Sec. II. If NSn-4-branes are condensed in the bulk, the link fields are 
confined. This is what happens in the confinement phase and the non-holographic critical phase 
discussed in the previous section. The bulk physics alone can not distinguish the confinement 
phase and the non-holographic critical phase. 

If the link fields are condensed, the two-form gauge field acquires a mass due to the Higgs 

nn 

mechanism ||43L I46Q. Note that loop fields with arbitrarily large size acquires expectation values 
in the Higgs phase because loop fields are just products of link fields. This corresponds to the 
deconfinement phase of the boundary matrix model. 
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If the link field is gapped and the two-form gauge coupling is small, the theory can be in the 
Coulomb phase. In this phase, closed strings are deconfined and the two-form gauge field arises 
as a light mode in the bulk. This corresponds to the holographic critical phase. Note that the loop 
fields can have finite expectation values even though link fields are gapped in this phase. 

E. Holographic critical phase II : deconfined open string 

The mean field description discussed in the previous section allows one to understand a yet 
another new phase of the matrix model. To see this, we first note that the decomposition in Eq. 
(l72l) has a U(l) gauge redundancy, 

£m« -> e l ^'-^£ M (i), (75) 

where ji is a U(l) phase defined on each site on the bulk space. Because of this U(l) redundancy, 
the link field can not have a quadratic hopping term. This is similar to the U(l) gauge redundancy 
present in the slave-particle theory discussed in Sec. II. One can decouple the quartic term for the 
link fields using the Hubbard-Stratonovich transformation. The resulting action should include a 
dynamical compact U(l) gauge field, 

S' bulk = -* £ [eM(^ + ^)eM(^)e- l ^ W -^ (i+M)+A - W1 +c.c.] 

i,M^N 

~t J2 [G(*)6f(^)e- i[i3Miv(i) ~ Aiv(i+M)+AM(i+Ar)] +c.c] 
-t £ [^(^ + M)af(0e" 4[i?Miv(l)+AM(i+Af)+Ajv(l)1 +c.c.] 

+ £^(KmW! 2 )-^- cos(A l B mn + A m B nl + A n B lm ). (76) 
iM 9kr cubes 

Under the gauge transformation in Eq. d75l) . the U(l) gauge field transforms as usual, A M {i) — > 
Am(i) + li+M — li, and the two-form gauge field is invariant. Eq. (|75T) implies that the link field 
£m(«) carries U(l) charge +1 on one end at site i + M and charge —1 on the other end at site i. 
Under the two-form gauge transformation in Eq. (1731) , the U(l) gauge field transforms as 

A M (i)^A M (i) + 9 M (i). (77) 

The first term in Eq. (1761 ) describes a link parallel to the direction M hops along the direction N 
as is shown in Fig. H3](a), and the second and third terms describe hoppings where the direction 
of a link changes as in Fig. H3](b) and (c). 
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i+N i+M+N i+N i+M+N i+N i+M+N 
> < > 



(a) 



i+M 



(b) 



i+M 



(c) 



i+M 



FIG. 13: The figures in (a), (b) and (c) show a 'particle' denoted as a thick line defined on the link 
< i,i + M > hops to the link < i + N,i + M + N >, < i,i + N > and < % + M + N, i + M >, 
respectively. The dashed lines represent the paths along which the end points of the link follow to which 
the U(l) gauge field is electrically coupled. 



Although the bare coupling for the U(l) gauge field is infinite, the kinetic term will be generated 
once high energy fluctuations of the link fields are integrated out, renormalizing the gauge coupling 
to a finite value. Suppose we are in the holographic critical phase where the link fields are gapped 
and the two-form gauge field is in the deconfinement phase. The U(l) gauge field can be either 
in the confinement phase or in the deconfinement phase. If the U(l) gauge field is confining, all 
open links are joined with each other to form closed loops because there is a linearly confining 
force between open ends. In this phase, only closed strings are allowed. This is the holographic 
critical phase I discussed in Sec. VII.C.3. On the other hand, if the U(l) gauge field is in the 
deconfinement phase, closed strings can get fractionalized into open strings, and open strings arise 
as deconfined excitations. Closed strings can still exist as a bound state of open strings. In this 
phase, there is an emergent U(l) gauge field in addition to the two-form gauge field. The two- 
form gauge field is coupled with the world-sheet of strings and the U(l) gauge field is coupled 
with boundaries of open strings. Here the gapless U(l) gauge field is protected by a quantum 
order associated with the emergent Bianchi identity dF = for the U(l) gauge field in the bulk. 
Open strings are fractionalized excitations of closed strings. 

The way open strings and the U(l) gauge field arise as collective excitations of closed string 
fields is very similar to the way slave-particles emerge as fractionalized excitations along with the 
emergent U(l) gauge boson in the slave-particle theory discussed in Sec. II. One can go one step 
further to obtain a different gauge group for open strings. For this we introduce a larger gauge 
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(a) 



(b) 



FIG. 14: (a) In the confinement phase for the one-form gauge field in the bulk, link fields are always joined 
to form closed loops, and there exist only closed strings as finite energy states, (b) In the deconfinement 
phase, closed strings can be fractionalized into open strings. The end points of open strings carry gauge 
charge for the emergent gauge field in the bulk. 



redundancy in Eq. (|72l) . 



>c 



» = tr 



n 

(i,i+tM)eC 



(78) 



where EJi, z) 



(i + fi, z) is a N x N complex matrix field defined on the link < i,i + jj, > 



in the bulk. This decomposition has the U (N) gauge redundancy, 



l^z) ->• V^(z)E^i,z)V i+ „(z). 



(79) 



Therefore the link fields now have to be coupled with a dynamical U (N) gauge field in the bulk. 



bulk 



-t 



E 



tr 



(i + N)Ui 



-iB M N(i) 



+ C.C. 



-t Y tr [^ N (i)E M (i)Ui + M,i+M+NU i+ M+N,i+Ne Wmn ® + C.C. 



-t tr[E N (i + M)E M (i)U i+M+N>i+N U i+N , i e- iBM ^ + c.c] 



+ Y V (\ E M(i) 



LM 



| 2 )-^- cos{A L B MN + A M B NL + A N B LM ). (80) 
9kr cubes 

Here U iti+M = e iTOi «®, and A a M (i) is U(N) gauge field defined on the link < i,i + M >. It 
is noted that the U (N) gauge field that emerges in the bulk is different from the original U (N) 
gauge field of the boundary matrix model. In this description, end points of the link fields carry 
fundamental and anti-fundamental charges for the dynamical U(N) gauge field. If the U(N) 
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gauge field in the bulk is in the deconfinement phase, open strings with the U (N) Chan-Paton 
factor emerge as collective excitations of the closed string fields. 

Although one can choose a description with an arbitrary gauge redundancy, the gauge group 
is determined by dynamics in the end. Holographic theories with different gauge groups for open 
strings in the bulk describe different states of the matrix model. In the T-dual description, configu- 
rations with background gauge fields describe D-branes. It is interesting to note that D-branes can 
emerge as non-perturbative excitations in the closed string field theory. 

VIII. DISCUSSION 



TABLE II: 





closed string 


open string 


NS-brane 


bulk excitations 


Confinement phase 


confined 


confined 


condensed 


X 


Non-holographic 
critical phase 


confined 


confined 


condensed 


X 


Holographic 
critical phase I 


deconfined 


confined 


gapped 


closed string (Bmn), 
NS-brane 


Holographic 
critical phase II 


deconfined 


deconfined 


gapped 


closed & open string (Bmn, Am), 
NS-brane, D-brane 


Deconfinement (IR free) phase 


condensed 


condensed 


confined 


non-local string 



In summary, we showed that a D-dimensional gauged matrix model can be mapped into a 
closed string field theory in (D + 1) -dimensional space. The string field in the bulk is coupled 
with a compact two-form gauge field which is also a part of the string field. Holographic states 
with deconfined string in the bulk are stable only when topological defects for the two-form gauge 
field are suppressed in the bulk, which is likely to be realized for a sufficiently large N. The holo- 
graphic states are in different universality classes from the non-holographic states where strings 
are confined in the bulk due to condensed topological defects. We also discussed a holographic 
critical state where closed strings get fractionalized into open strings. In this state, there are both 
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closed and open strings along with the two-form and one-form gauge fields in the bulk. The non- 
trivial quantum order present in the holographic phases is responsible for the existence of operators 
whose scaling dimensions are protected, which otherwise would have received a large quantum 
correction at strong coupling. The possible phases of the matrix model are summarized in Table. 

n. 

Although many structures on the holographic description have been learned from general con- 
siderations, it is desirable to obtain explicit solutions to the saddle point equation. In principle, 
one has to solve a set of infinitely many coupled differential equations. In the future, it will be 
interesting to simplify these equations by focusing on light modes. Finally, we close with discus- 
sions on a speculative phase diagram of the matrix model, a world-sheet description of deconfined 
strings in the holographic phases, and a continuum limit. 

A. A schematic phase diagram 

It may be difficult to find a specific microscopic model which realizes each phase discussed 
in this paper. However, one may still guess a possible phase diagram. For D < 4, it is believed 
that the present matrix model is always in the confinement phase. In these low dimensional cases, 
one may have to introduce more degrees of freedom (fermions or fundamental matters) to stabi- 
lize critical phases. Here we focus on the pure bosonic matrix model in D > 4. In Fig. \T5\ we 
show a speculative phase diagram. In the strong coupling limit, the matrix model is in the con- 
finement phase. As the gauge coupling is weakened, there is a phase transition into the IR free 
deconfinement phase. If the phase transition is continuous, there can be two different universality 
classes for the critical points. For small N, iVS£>_4-branes are condensed and strings are confined 
in the bulk. At this non-holographic critical point, the scaling dimension of Wilson loop operators 
generally receive a large quantum correction. For iV greater than a critical value, iVS^-4-branes 
are suppressed, and strings are deconfined in the bulk. The quantum order protects the two-form 
gauge field from acquiring a large mass in the large iV limit, which in turn protects the scaling 
dimension of the phase fluctuations of Wilson loop operators even at large 't Hooft couplings. 

As was discussed in Sec. VII. E, there exist two different kinds of holographic critical points. In 
the first case, there are only closed string excitations in the bulk. In the second case, there are both 
closed and open string excitations along with the emergent gauge field and the two-form gauge 
field. We expect that it is easier to stabilize the state with both closed and open strings for iV's that 
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FIG. 15: A proposed phase diagram for the pure bosonic gauged matrix model in D > 4. Here A represents 
a set of 't Hooft couplings associated with various Wilson loops in the matrix model. Moving along the 
direction of A in this figure may mean tuning more than one microscopic parameters at the same time. The 
critical points, in general, represent multi-critical points. Therefore the shape of the phase boundary in the 
figure should not be taken seriously. The actual phase boundary is not likely to be a straight line in the 
multi-dimensional space of microscopic couplings. In the strong coupling limit (large A), the theory is in 
the confinement phase with the exponentially decaying Wilson loop correlation function. As the coupling 
becomes smaller, the model goes through a phase transition to the IR free deconfinement phase. For small 
N, A r S*£)_4-branes remain condensed at the critical point. Strings are confined, and the two-form gauge 
field is massive in the bulk. Critical correlations are mediated by the D-dimensional fluctuations of strings 
near the UV boundary. Scaling dimensions generally receive a large quantum correction at strong couplings. 
This is the non-holographic critical point. For large N, A5D_4-brane is gapped and the two-form gauge 
field remains light even at strong coupling. Low energy physics is governed by deconfined strings that 
propagate deep inside the bulk. In this holographic critical point, critical correlations of the matrix model is 
mediated by (D + 1) -dimensional fluctuations of strings. The scaling dimension of the operator associated 
with the two-form gauge field is protected from acquiring a large quantum correction due to the quantum 
order associated with the dynamical suppression of A^S'D-4-branes. The holographic critical point can be 
either in the state with only closed strings, or in the state with both closed and open strings. 



are large but not too large. In the large N limit, closed strings are free, and there is no dynamical 
reason why they decay into open strings. For A^ large enough to suppress A r S , £)_ 4 -brane, but still 
small enough to support strong interactions between closed strings, closed strings may decay into 
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open strings. We believe that further studies are needed to understand this phenomenon more 
systematically. 

It is of note that the structure of the proposed phase diagram is reminiscent of known examples 
where systems flow into novel universality classes at interacting critical points. For example, in the 
two-dimensional Z N clock model with N greater than a critical value, the critical point between 
the disordered phase and the ordered phase has an emergent £7(1) symmetry [48]. More recently, 
it has been proposed that the critical point between an antiferromagnetic state and a valence bond 
state in 2+1 dimensions can possess a non-trivial quantum order which supports an emergent gauge 
boson and fractionalized excitations [49]. 

B. World sheet description of deconfined string 

In order to make a contact with the traditional first quantization formulation of string theory, it 
will be useful to have a world-sheet description of deconfined strings in the holographic phases. 
Here we focus on closed string. The generalization to open string is straightforward. Note that the 
hopping integral from loop G\ to C\ + C2 is determined by the loop field 4>c 2 = \<i>c 2 \e lbc ' 1 which 
is complex. The amplitude \<pc 2 | determines the strength of the hopping, and defines the notion of 
'distance' between the two loops. The distance between two loops, in turn, defines the metric of 
the space in which loops are defined. In this sense, condensates of loop fields determine the metric 
of the space in which closed loops propagate. The U(l) phase b C2 corresponds to the background 
two-form field to which closed strings are electrically coupled. 



FIG. 16: The world sheet action has two contributions : A z is the area of the surface which is parallel to z, 
and Ac is the contribution from hoppings mediated by condensates of loop fields. 
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This can be made more intuitive if we use a world-sheet representation. Let us consider the the 
quadratic Hamiltonian that includes the tension and the hopping terms, 

H (z) = LcXcXc - ^^c a ( z )Xa l Xc 1 +c a + 24>c 2 (z)Xc 1+ c 2 Xc^ ■ (81) 

Here we suppressed the form factors F i j(Ci,C 2 ) and Gij(Ci, C 2 ). For this discussion, we ig- 
nore a possible deformation of the path integral of the string fields. We consider the single loop 
propagator given by 

g(C 2 ,z 2 ;C 1 ,z 1 ) = < C 2 \e~ ^ H ^ dz \C x >, (82) 

where \C >= Xc\® > is a single string state. We can 're-discretize' the imaginary time z into 
small steps with size e to write 

g(C 2 ,z 2 ;C 1 ,z 1 ) = £ e" 5 -^ 5 ^, (83) 

all world sheets 

where 

S^ s = A z - ^m(2e|0 c ^)|/M 2 ), (84) 

i 

where A z is the area of the world sheet whose faces are tangential to the z direction, and the second 
term in S^ s includes contributions 2e ^^ Zl ^ for each hopping mediated by loop field §q at ti me 
Zi. The second term is associated with the parts of the world sheet that are perpendicular to the z 
direction. This is illustrated in Fig. [TH One can view S l sw as the Nambu-Goto action provided 
that the area associated with a loop C at time z is taken to be A c = — \n(2e\^c( z i)\ I M 2 ) . The 
areas associated with loops in turn would determine the spatial metric of the space. The imaginary 
part of the action, 

S 2 WS = [ B (85) 

./world sheet 

is simply the Berry phase associated with the phase of the background loop fields in the temporal 
gauge. Since both the phase and amplitude modes of loop fields are dynamical, not only the 
two-form gauge field but also a metric field should arise as a dynamical degree of freedom. 

A typical configuration of vacuum fluctuation in the full interacting string theory are shown in 
Fig. [TTJ For every vertex where closed loops join or split, there is a factor of 1 /N. For a large N, 
the theory describes weakly interacting strings propagating in the time-dependent background. 
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FIG. 17: A snap shot of a vacuum fluctuation in the holographic critical phase with deconfined closed 
string. 



C. Continuum limit 

In the holographic phase, there are three important length scales. The first scale is associated 
with the tension ^ns D -4 °f NSD-A-brane, Ins D - 4 ~ ^-/(^NSn-J 3 - The second scale is the 
'string scale' l s which corresponds to the typical size of closed string excitation in the bulk. The 
third scale L is the scale over which loop fields change appreciably in the bulk. Roughly, the last 
one determines the 'curvature' of the bulk space in which strings propagate. One can take the 
continuum limit by tuning the system such that all these length scales are fixed in the a — > limit. 
If these scales satisfy L » l s » Ins d ^ 4 , strings propagate in a weakly curved background. 
It is expected that a continuum description of string theory emerges in this limit. Since graviton 
has the same mass as the two-form gauge field in the continuum limit, graviton may emerge as a 
massless mode along with the two-form gauge field in the holographic phase. It will be interesting 
to see how the resulting theory in the continuum limit compares with the existing formulation of 



the closed string field theory[50] 
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X. APPENDIX A 



We prove the identity 

I = J d<p*d<P e-^- A * } /(0*) 

= tt/(A*) (86) 

for any analytic function f(x). We define a new variable £ = — A to write 

I = J d£*d£ e~^ 2 - A ^f(C + A*) 

= EEti / d ^ e-W(-AC) m f in \A*)r, (87) 

where we Taylor expanded e A ^* and /(£* + A*). Here only m = n = component survives in the 
angular integration of £ = \£,\e ld and we have 

/ = nf(A*). (88) 

XI. APPENDIX B 

Here we show that Eq. (1391 with Eqs. (|40l ) - (|43l is equivalent to Eq. (|34l ). The imaginary time 
j3 is divided into R pieces, 

, R-l 



Z = ff[d*<§*d*® <^ / |e-^|$^ 1) >e-*c- 1) ^%$g ? - 1 )|e^|$ 
J i=i 



(R-l) la -eH l(h (R-2) > x 



e-^c- 2) ^c R - 2) ... < $g)| e -*|$g) > e"*^ 15 < $gV^ >, (89) 
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where we use the identity 



For (|40l ) and (|42l) . the transition amplitude becomes 



J 1=0 



where 



R 



We take i? — >• oo limit and equate Eq. (|92b with Eq. (1341) to identify 



(90) 



(91) 



(92) 



(93) 



XII. APPENDIX C 















<t>n 






4>n 







FIG. 18: 



In this section, we solve the loop equations perturbative in 1/M to the lowest non-trivial order 
for small loops. Let us first consider the saddle point solution in the large M limit, which corre- 
sponds to the deep inside the confinement phase of the gauge theory. For M = oo, we have simple 
solution, 



<f>c{z) = <M0K 



-Lqz 



<j>* c {z) = #,(0)e 



Lqz 



(94) 

with the boundary conditions in Eqs. (|651) and (|66l) . Especially, the IR boundary condition implies 
4>q(oo) = 0, which then implies <p* c {z) = 0. This makes sense physically because the source for 
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Wilson loops decreases exponentially with z, and the expectation values of Wilson loop operators 
vanish in the strong coupling limit of the gauge theory. For a large but finite M, the loop equations 
can be solved perturbatively in 1/M. 

To illustrate the idea, we solve the loop equation for a simple model which include only the 
single trace operator for the unit plaquette, V = J n <p* D . In this model, we solve 0* to the lowest 
order in 1/M. Keeping the first four shortest loops that contribute to the evolution of <p* c of unit 
plaquette, we write the saddle point equations, 

2 

<9 2 0* = 40* - j^(4>* u + 4>n + 0c + 0c)0s> 
d z 6* = 60*,--4r(0T + 



1*2 



= ^l-Jp^l +£-). 



«9 2 0_ = 20* - — . (95) 



2 

M 2 ' 

The symbols are defined in Fig. [18] We consider solutions with the translational and discrete 
rotational symmetries of the lattice, and suppress the site index for the loop fields. The equations 
for 0p, 0*_ and D are similar to the one for 0*j. Similarly, the equation for 0* (0|) is similar to the 
one for 0^ (0-). From Eqs. (|32~1) and (|66l) . we read the boundary conditions, 

,*/ x 0S(OO) , . 1 

Using the zero-th order solution for Q = J Q e _4z , we obtain 



M 8 ' 



•□(*) = T^e" 42 . (97) 



This implies that the expectation value of the Wilson loop operator for the unit plaquette become 
Jg/M 8 as expected in the large M limit. The difference from the case with M = oo is that <f)* c is 
now non-zero at z = because Wilson loops have non-zero expectation values. As expected, 0c 
and (j)* c satisfy Eq. (1671 to the leading order in 1/M at all z. 
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